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With a Geant4 software package based on the Monte Carlo method, a multi-cell 47t detection system is de- 
signed, which consists of 40 Gadolinium-loaded liquid scintillation detectors. These detectors, associated with 
a fission chamber in its geometrical center, constitute a platform. This platform is mainly used for the measure- 
ment of a fissionable nucleus (n, 2n) reaction cross section. In order to properly determine the experimental 
set-up, we carry out a systematic numerical simulation using our model which is established by the Geant4 soft- 
ware package. This work provides rich and valuable reference data for experiments on the fissionable nucleus 


(n, 2n) cross section measurement in the future. 
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I. INTRODUCTION 


(n, 2n) nuclear reaction plays the role of neutron multi- 
plication in the reactor, so an accurate measurement of the 
fissionable nucleus (n, 2n) cross section has an important 
significance in the utilization of nuclear energy and nuclear 
technology [1]. People have employed a variety of methods 
to carry out this work. Activation technique is one of the 
commonly used methods, but it requires the residual nucleus 
appropriate half-life. In addition, this method has a higher 
requirement for the purity of the samples and the intensity 
of incident neutrons [2, 3]. Partial y-ray method is a tech- 
nique that combines with the theoretical model loading to a 
nuclear reaction. But it depends too much on the theoretical 
model, which leads to great uncertainty to the experimental 
results [2, 3]. The large Gadolinium-loaded liquid scintillator 
was introduced by Frehaut et al. in 1976 for its 80% neutron 
detection efficiency [4, 5]. The principle of this direct mea- 
surement method is shown in Fig. 1, two hemispherical con- 
tainers contain liquid scintillator and a dozen photomultiplier 
tubes surround the containers to collect the signal from the 
Gd-liquid scintillator. After being collimated, the monoener- 
getic neutrons irradiate the samples located in the geometric 
center of the system. The neutron coming from the nuclear 
reaction is moderated within a few nanoseconds in the liquid 
scintillator, and generates a neutron pulse (usually called fast 
signal). Since Gd isotopes mainly !??Gd and !5"Gd have a 
high thermal neutron capture cross section, there is a great 
probability that the neutron will be captured by the Gd. Gd 
isotopes in an excited state will emit several cascade y-rays 
with a total energy of about 8 MeV, and the pulse caused by 
these captured y-rays is called the slow signal. Because of 
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Fig. 1. Experimental set-up. 


the high energy of this pulse, we can easily discriminate it 
from common background events such as y-ray and scatter- 
ing neutron events, whose energy is less than 3 MeV. And 
such a pulse with an energy of 8 MeV indicates that a neu- 
tron is probably produced in the nuclear reaction. Thus, we 
can exactly determine the number of neutrons produced in a 
nuclear reaction. 


This method does not rely on the theoretical model, and has 
no requirement for the half-life of the residual nucleus [4—7]. 
However, the fission neutrons will bring great deviation to the 
measuring results of the fissionable nucleus (n, 2n) cross sec- 
tions. In addition, because the detection system has a large 
sensitive volume, the y-rays in the background have a signif- 
icant additive effect, resulting in a high y background. And 
each y-ray penetrated in the detection system can make all 
the photomultipliers respond, so the total counting rate of the 
detection system is very high, which results in a high dead 
time for the detection system. 

Considering the above disadvantages of the large Gd- 
loaded liquid scintillator technique, we improved this tech- 
nique by a novel design of one 47t detection system, which 
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Fig. 2. (Color online) Profile of the detection system. 
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Fig. 3. (Color online) Transverse cross-section of the detection sys- 
tem. 


consists of 40 Gd-loaded liquid scintillatior sub-detectors 
(Figs. 2 and 3). In this experiment, the multi-cell structure 
can effectively reduce the counting rates of a single detector, 
so as to decrease the dead time. We will also see from Sec. 
III A of this paper that this technique has certain (n, y) dis- 
criminations, which can effectively reduce the influence of y 
background on the experimental results. To measure the fis- 
sionable nucleus (n, 2n) cross section, samples are placed in 
the fast fission chamber, which is located at the geometrical 
center of the detection system. The pulse output from the 
fission chamber will be used to exclude the fission neutron. 
Thus, the influence of the fission events on the experimental 
results will be greatly reduced. This platform can not only be 
used to measure the (n, 2n) and (n, 3n) cross sections, but also 
to measure the neutron properties such as the angular distri- 
bution and angular correlation. Thus, the results provided in 
this work are very important for neutron experiments in the 
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future. 


Il. STRUCTURES AND BASIC PERFORMANCE OF THE 
DETECTION SYSTEM 


A. Structure of the detection system 


In order to cover the 47 solid angle, every two adjacent de- 
tection units in the detection system should realize the seam- 
less connection. At the same time, each detection unit should 
have an approximately equal solid angle to the geometrical 
center of the system to maintain the same detection efficiency, 
According to Ref. [8, 9], a detection system containing 40 de- 
tection units is designed, in which each unit is numbered as 
shown in Fig. 2. The detection units used in this system have 
two kinds of shapes whose transverse sections are pentago- 
nal and hexagonal, respectively. In the center of the detection 
system a cavity is left (Fig. 3), which is used to hold the fast 
fission chamber. Furthermore, in order to design the neutron 
beam tube, two detection units are extracted from the detec- 
tion system (Fig. 3). Each detection unit is a container which 
has an aluminum shell with a thickness of about 3 mm. These 
detection units are all filled with the common liquid scintil- 
lator used in the neutron experiment. The liquid scintillator 
is mixed into a Gd isotope by a mass percentage concentra- 
tion of 0.596. The ratio of carbon to hydrogen is approxi- 
mately 3 : 4 in the liquid scintillator and the density is about 
0.87 g/cm?. 


B. Performance of the detection system 
1. Response of the detection system to the neutron source ?? Cy 


The spontaneous fission neutron source, 2520 f> is often 
used for the efficiency calibration in the neutron detection ex- 
periment. Therefore, we simulate this kind of neutron source, 
put it in the geometric center of the detection system, and ana- 
lyze the response of the detection system to it with our Geant4 
model. Partial y-rays are emitted from excited Gd isotopes 
after they capture a neutron, we show the energy spectrum 
that these y-rays deposit in detection unit and detection sys- 
tem, respectively (Fig. 4). Since y-rays primarily lose their 
energy in the form of the Compton Effect in a liquid scintilla- 
tor [10], the y-ray energy spectrum of a single detection unit 
is the Compton plateau whose edge is approximately at the 
position of 8 MeV. The detection system has a large sensitive 
volume, and the energy of these cascade y-rays can be fully 
deposited in it after several Compton scatterings, thus a full- 
energy peak is formed at the position of 8 MeV. In addition, 
we also show the neutron efficiency curve in Fig. 5 (In Figs. 5 
and 6, ‘Th’ means the threshold value of a single detection 
unit). The neutron efficiency is defined as 
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Fig. 4. (Color online) y-rays spectrum of the detection system and a 
detection unit. 
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Fig. 5. (Color online) The curve of the neutron detection efficiency 
versus the thickness of the detection unit. 


where c is the neutron efficiency, Ny is the total counts in 
the y-ray spectrum caused by neutron events whilst the JV;, 
is the total number of the neutron events. In the experiment, 
the detector threshold is usually far below 1.0 MeV. We can 
see from Fig. 5 that the efficiency of a single neutron is over 
8096, except for the curve at the bottom, and they increase 
rapidly with the increasing thickness of the detector. When 
the thickness is over 60cm, the efficiency curve gradually 
flattens out. It is indicated from Fig. 6 that the variation of 
the cavity radius has limited influence on the coincidence ef- 
ficiency of two neutrons, except when the detector threshold 
is 2.0MeV. To minimize the volume of the whole detection 
system, which also has a high neutron efficiency, the thick- 
ness of the detector is preliminarily determined to be 60 cm. 
Furthermore, considering a fission chamber will be placed in 
the center of the detection system, we chose the cavity radius 
to be 30cm. 
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Fig. 6. (Color online) The curve of coincidence efficiency about two 
neutrons versus the inner-radius of the detection system. 
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Fig. 7. (Color online) Statistical results about the distribution of the 
energy deposition in detection units. 


2. Distribution of y-ray energy deposition in the detection system 


The neutron source, ?9?C +, which is placed in geometric 
center of the detection system emits neutrons in the direction 
of the detection unit 36 (Fig. 2). If the neutron is captured 
by the Gd isotope in unit 36, the deposited energy of the 
capture y-rays in each detection unit will be recorded. Ac- 
cording to the simulation, we give out the distribution of the 
energy deposition in the detection system, which is shown in 
Fig. 7. We can see in Fig. 7 that the y-rays deposit most of 
their energy in several adjacent units such as 36, 20, 21, 26, 
11, and 27. These adjacent detection units are called partial 
detection units, where about 60% of the captured y-rays en- 
ergy is deposited. This characteristic can be understood as 
follows. Although the total energy of the captured y-rays is 
close to 8 MeV, they commonly include several y-rays whose 
energy is about 1~2 MeV. Because their penetrating power is 
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very limited, most of their energy is deposited in several adja- 
cent detection units. This indicates that partial detection units 
(composed of 6~7 adjacent detection units) will respond to 
the event fired by a neutron. However, each y-ray coming 
from the background is commonly low-energy, so it will de- 
posit most of its energy in one or two detection units, which 
is obviously different from the neutron event. Thus, the num- 
ber of fired sub-detection units can be used as a judgment of 
a neutron event in the future, so as to reduce the y-rays back- 
ground coming from the environment. 


Ill. NUMERICAL SIMULATIONS ABOUT THE 
EXPERIMENT 


A. Experimental set up in simulation 


The fast fission chamber used in the detection system is a 
multi-target chamber, which can effectively reduce the self- 
absorption of the fission fragment without influencing the 
(n, 2n) reaction rate [11]. The (n, 2n) reaction threshold of 
the aluminum nuclei is relatively higher than other common 
metallic elements. Even if the (n, 2n) channel is induced, 
its cross section is also small. So, duralumin is chosen as the 
structure material of the fast fission chamber [11, 12]. The to- 
tal mass thickness of the structure material is 2.7 g/cm?. The 
fission chamber is placed in the geometric center of detec- 
tion system and its axis coincides with the axis of the neutron 
beam tube. ?35U is plated on both sides of the electrodes 
in the chamber, and the total mass thickness of the sample 
is about 100 mg/cm?. The position of the monoenergetic neu- 
tron source (14 MeV) is 3.5 m away from the geometric center 
of the detection system. With these conditions, the incident 
neutron may induce many kinds of nuclear reactions. 


B. Characteristic of time distribution of the different events 


In order to properly set up the pulse period of the incident 
neutron, the width of the detector gate and the time relation- 
ship between the gate signal and the neutron flux, we simulate 
the time distribution of different events which induce the re- 
sponse of the liquid scintillator. These events mainly include 
the fast signal (elastic scattering neutron, inelastic scattering 
neutron, fissional neutron, and partial y-rays which are also 
produced by the samples) and slow signal (captured y-rays 
which are emitted from excited Gd isotopes). In our sim- 
ulation, we chose the time that the neutron source emits a 
neutron as the zero time. The time distribution of these dif- 
ferent events is shown in Fig. 8. As shown in Fig. 8, what first 
causes the detector to respond is the prompt y-ray, and then 
the neutron events. The time distribution of these events is in 
the nanosecond order of magnitude; while Figs. 8 and 9 show 
that the time distribution of the slow signal is in microsecond 
magnitude. Our simulation also shows that 99.5% of the fast 
signals distribute within 1 us. And only about 0.2% of the 
slow signals distribute in this interval (blue curve in Fig. 8). 
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Furthermore, 95% of the slow signals distribute within 30 us 
and more than 99.7% distribute within 50 us. 
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Fig. 8. (Color online) Time distribution of different events including 
fast and slow signals. 


It is noted that the time distribution of slow signal is ap- 
proximately equal to that of neutron capture by Gd. And this 
slow signal is what we are actually interested in, as it is actu- 
ally used for neutron counting. To avoid the influence of the 
fast signal, the detector gate should be opened after about 1 us 
of the neutron flux arrival. The time width of the detector gate 
signal should be set to about 30 us, during which the neutron 
is counted. The time between two consecutive neutron fluxes 
is chosen to be 60 us, which is greater than the lifetime of neu- 
trons in the scintillator. In Ref. [4—6], the above parameters 
used by Frehaut et al. are 1.2 us, 30 us and 60 us, respectively. 
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Fig. 9. Time distribution of slow signal, which reflect the time dis- 
tribution of neutron capture by Gadolinium. 
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TABLE 1. The numbers of different 2n events which induce the response of the detection system in our simulation 


235U(n, 2n) (n, f) 27 Al(n, 2n) MN(n, 2n) Accidental Coincidence Total Counts 
Number of Events 9411 750 833 2142 643 13779 
Counting Uncertainty 27.36 60.15 117.38 


C. Estimation about the accuracy of the results 


With a total number of 1.0 x 10? neutrons irradiating the 
samples in the fission chamber, we carry out a systematic sim- 
ulation of the experiment. Table 1 lists the kinds of the 2n 
events that the detection system records. These 2n events can 
be ascribed to the following reasons, the fission reaction of 
uranium, the (n, 2n) reaction of uranium, aluminum, and the 
accidental coincidence. The (n, 2n) reaction of nitrogen (the 
main element in the air) should also be considered. If the ef- 
ficiency of the chamber is about 90%, only 10% of the 2n 
events (this number is about 750 in our simulation) from the 
fission reaction can not be excluded from the experimental re- 
sults. That is to say, the influence of the fission neutron on the 
experimental results can be reduced by one order of magni- 
tude. When calculating 2n accidental coincidence events, the 
coincidence resolving time is equal to the width of the detec- 
tor gate signal (30 us); the total intensity of the neutron source 
in the 47 direction is 105/s (The 10°/s~10°%/s intensity can be 
achieved on the accelerator at the China Institute of Atomic 
Energy). 

The 2n background events caused by the structure materials 
(aluminum) of the chamber and nitrogen in the air can be de- 
ducted through the following method. Under the same condi- 
tions, we can determine these 2n background events with the 
neutron beam irradiating the fission chamber without sam- 
ples. Furthermore, simulation also shows that only 1096 of 
the accidental coincidences are related to the samples [13], 
namely, 9096 of accidental coincidences can also be deducted 
through the above method. Therefore, the counting uncer- 
tainty of 2n events caused by aluminum, nitrogen, and acci- 
dental coincidence can be calculated together. And the uncer- 
tainty of net count can be estimated 


a= joi +0} toj = 135. (2) 


Here, c; = 117.38, op = 60.15, of = 27.36 which are 
listed in Table 1. Relative uncertainty is 


on 


= 1.43%, (3) 


Or = 
net 
where Nnet = 9411 is what we are really interested in. 
Of course, a few factors such as the accidental coincidence 
with cosmic rays and the energy resolution of the y-ray spec- 
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trum are not considered at present. There will be some devi- 
ation in calculating the content of the sample; and the change 
of the neutron flux can also lead to deviation. So, the mea- 
surement of the (n, 2n) cross section will be influenced by 
many factors, except for those which we have considered. 
Correspondingly, the value of the experimental result uncer- 
tainty will be larger than that we estimate at present. How- 
ever, the simulation results in our study are still important to 
the design and optimization of the experiment conditions in 
the future. 


IV. CONCLUSION 


According to the experiment setup of Frehaut et al., we 
designed a novel 47t multi-cell Gd-loaded liquid scintillator 
detection system for the measurement of (n, 2n) nuclear re- 
action cross sections. This detection system makes up for the 
weakness of the original one used by Frehaut et al. According 
to systematic simulation of the experiment, we preliminarily 
determined the thickness of the detector, the inner-radius of 
the detection system, the pulse period of the neutron source, 
the time width of the detector gate, and the time relationship 
between gate signal and neutron flux. Finally, we roughly 
estimated the accuracy of the experimental result under our 
experiment setup. 


Certainly, there is much work that needs to be further car- 
ried out. As for the y-rays spectrum, we only considered the 
energy loss process in the detector in our simulation, without 
considering the collection process of the fluorescence pho- 
ton. Considering this process, the y-rays spectrum will have 
a lower energy resolution, and the energy resolution has an 
important influence on the efficiency. Therefore, the photo- 
gathering process should be considered in our Geant4 code in 
the future, so that the simulation will be more consistent with 
the experimental conditions. In addition, in order to reduce 
the influence of the large angle scattering neutrons, cosmic 
rays and y-rays, coming from the environment, we will put 
forth much effort to design a neutron beam tube and a shield- 
ing layer around the detection system. It mainly includes the 
selection of shielding materials, the calculation at the thick- 
ness of the shielding layer, and the diameter of the neutron 
tube. 
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